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ABSTRACT 

We report the discovery of a massive ultra-compact quiescent galaxy that has been strongly-lensed 
into multiple images by a foreground galaxy at z = 0.960. This system was serendipitously discovered 
as a set of extremely K s -bright high-rcdshift galaxies with red J - K s colors using new data from the 
UltraVISTA YJHK S near-infrared survey. The system was also previously identified as an optically- 
faint lens/source system using the COSMOS ACS imaging by Faure et al. (2008, 2011). Photometric 
rcdshifts for the three brightest images of the source galaxy determined from twenty-seven band 
photometry place the source at z = 2.4 ± 0.1. We provide an updated lens model for the system 
which is a good fit to the positions and morphologies of the galaxies in the ACS image. The lens model 
implies that the magnification of the three brightest images is a factor of 4 - 5. We use the lens model, 
combined with the K s -band image to constrain the size and Sersic profile of the galaxy. The best-fit 
model is an ultra-compact galaxy (R e = 0. 64io.il kpc, lensing-corrected), with a Sersic profile that 
is intermediate between a disk and bulge profile (n = 2.2^o 9), albeit with considerable uncertainties 
on the Sersic profile. We present aperture photometry for the source galaxy images which have been 
corrected for flux contamination from the central lens. The best-fit stellar population model is a 
massive galaxy (Log(M star /M ) = 10.8+0^, lensing-corrected) with an age of 1.0^;° Gyr, moderate 
dust extinction (A v = 0.81q'|), and a low specific star formation rate (Log(SSFR) < -11.0 yr _1 ). This 
is typical of massive "red-and-dead" galaxies at this redshift and confirms that this source is the first 
bona fide strongly-lensed massive ultra-compact quiescent galaxy to be discovered. We conclude with 
a discussion of the prospects of finding a larger sample of these galaxies. 

Subject headings: infrared: galaxies - galaxies: evolution - galaxies: high-redshift - galaxies: structure 
- galaxies: fundamental parameters 



1. INTRODUCTION 

The brightening and magnifying effects from strong 
gravitational lensing are a powerful tool that permit 
us to study high-redshift galaxies with better signal- 
to-noise (S/N) and higher spatial resolution than is 
normally possible with current instrumentation. In 
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the last decade significant successes in finding samples 
of lensed galaxies have come fro m surveys that have 
searched behind clusters (e.g., iBavliss et al.l 1201 lblj 
and luminous red galaxies (e . g.. | Bolton et al.l 120081 : 
iFaure et al.l 120081: ICoorav et aLl l2011h . These samples 
are now giving us insights into the evolution of typical 
galaxies at 1 < z < 5 which are mo st frequently low- 
mass blue star- forming galaxies (e.g.. | R i gbv et all 120111: 
iWuvts et al.1 l2012t ISharon et al.l l2012t iBrammer et al.l 
I2012D With the advent of the Herschel and South 
Pole Telescope sub-millimeter observatories we have 
also gained access to a population of strongly-lensed 
submillimeter galaxies which are readily identifi- 
able in wide-field surveys (e.g., iNegrello et aLl [2 010: 
I Vieira et al.ll2010|) or behin d clusters (e.g.. lEgami et all 
l2010t iGladders et all 1 2011) . These lensed submillimeter 
galaxies are providing intriguing results about how 
dusty and str ongly star forming galaxies evolve at high 
redshift (e.g., ISwinbank et~al1l2010l 1201 ll) . 

What has remained elusive are strongly-lensed 
examples of the proverbial "red-and-dead" massive 
galaxies at z ^2. These galaxies have generated 
significant interest in the last few years after it was 
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Fig. 1. — Postage stamp optical/NIR images of the lens/source system with a high-contrast stretch. The field-of-view in each image is 
12" X 10". The lensed galaxy is faint in the observed optical, and has an extremely red J - K s color. The separation between the multiple 
images is ~ 2". 

deblended from contamination from the central source. 
In § 5 we use this photometry to determine an accu- 
rate photometric redshift as well as stellar population 
parameters for the galaxy. We conclude in § 6 with a 
summary and a discussion and prospects of finding more 
strongly-lensed massive ultra-compact quiescent galaxies 
in future deep, wide-field NIR surveys. Throughout this 
paper we assume a fl\ = 0.7, £l m = 0.3, and Ho = 70 
km s _1 Mpc^ 1 cosmology. All magnitudes are in the 
AB system. 



iBruce et al.ll2012L and numerous others). Understanding 
how they grow so substantially in size without adding 
significant stellar mass is currently one of the major 
open questions in galaxy evolution theory. 

Finding gravitationally-lensed examples of these 
sources could prove to be valuable for understanding 
their evolution. The magnifying effects from lensing 
will allow us to potentially resolve the central regions of 
these compact sources and understand what their central 
stellar density profiles are. The brightening effects from 
lensing will also allow us to obtain higher S/N spec- 
troscopy of these systems. Currently, even the brightest 
of these systems are extremely faint and determining 
quantities such as velocity dispersions and dynamical 
masses requires substantial integrat ion times even on 
the world's largest telescopes (e.g., Ivan Dokkum et afl 



20091 ICappellari et al.l 120091: Ivan de Sande et all l201lT 
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Unfortunately, these galaxies are much more challeng- 
ing to detect as lensed sources and until now none have 
been securely identified. The challenge is that while they 
comprise ~ 50% of the population o f massive galaxies 
(LogM. strar /M^ > 11 0) at z ~ 2 (e.g.. lKriek et al.ll2008l : 
iBrammer et al.l 1201 ll ). they are still quite rare (~ 1 
every 10 arcminute 2 ), making it unlikely to find one in a 
favorable alignment with a foreground lensing structure. 
Furthermore, they are extremely faint in the observed 
optic al bands and also have very red J - K s NIR colors 
(e.g., iFranx et al.l 120031 : iKriek et al.ll2008ft making both 
deep and wide-field NIR imaging a requirement for their 
detection. 

In this paper we report the discovery of the first exam- 
ple of a massive ultra-compact quiescent galaxy that has 
been strongly lensed by a foreground galaxy. The galaxy 
was serendipitously discovered as a set of bright and 
extremely red high-redshift galaxies using new data from 
the 1.8 deg 2 UltraVISTA YJHK, near-infrared (NIR) 
survey (see iMcCracken et al.l [2013) . The lens/source 
system was previously identified as an optically faint 
quadruply-lens ed galaxy in the COSM OS field using 
ACS imaging (jFaure et all 120081 1201 ll) : however, with 
the addition of the UltraVISTA NIR data it is now clear 
that the source is a massive quiescent galaxy. 

This paper is organized as follows, in § 2 we discuss 
how the lens/source system was discovered. In § 3 
we present an updated lens model for the system and 
use this to constrain the size and Sersic profile of the 
source galaxy. In § 4 we present updated twenty-seven 
band photometry for the source galaxies that has been 



2. IDENTIFICATION OF THE LENS/SOURCE SYSTEM 

The lens/source system was serendipitously discov- 
ered using photometry from a K s -selected catalog of 
the COSMOS/UltraVISTA field. The catalog contains 
photometry in twenty-seven photometric bands includ- 
ing the publically-available seven optical broad band 
(u*g + r + i+z + BjVy) and twelve optical medium band 
(IA427 - IA827) i maging of the COSMOS field from 
iCapak et al.l ((2007). It also includes the public- release 
YJHK S NIR imagin g of the COSMOS fiel d from the 
UltraVISTA survey (jMcCracken et al.ll2012f ). as well as 
the public four-c hannel IRAC imagi ng is from the S- 
COSMOS survey (jSanders et al.|[2007l) . Object detection 
for the catalog was performed in the K s -band, which 
has a 5a depth in a 2" aperture of K s — 23.9 AB 
(jMcCracken et ai1 l2012). Ph otometry was measured us- 
ing the Sextractor package (jBertin fc Arnoutslll996D in 
2.1" apertures using PSF-matched images in all bands. 
Full details of the COSMOS/UltraVISTA catalog will be 
presented in a future paper (A. Muzzin et al., in prepa- 
ration). 

The system (R.A., 10:00:50.55, Decl, +02:49:01, 
J2000) was identified as a strong -lensing system via eye- 
examination of the ACS F81 4W (iKoekemoer et "al"1l2007ft 
and UltraVISTA K s -band (jMcCracken et al.ll2012f )"im- 
ages. This process was performed as a quality check of 
the catalog for a subsample of galaxies that were were 
bright (K s < 21.0 AB), red (J - K s > 1.5 AB) and had 
photometric redshifts (z p h oto ) > 2. In Figure 1 we plot 
postage stamp images of the system in the ACS F814W, 
and UltraVISTA J and K s bands with a high-contrast 
stretch. Sources A and C (see labels in Figure 1) were 
the sources in the catalog that matched the selection cri- 
teria. Sources B and D lie along the major axis of the 
lens galaxy and hence are more blended with the lens and 
were not identified as unique sources by SExtractor. The 
extraordinary brightness of sources A and C given their 
Zphoto, combined with the distinct cross-shape pattern of 
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four sources around the brighter central lens (similar to 
quasar strong lenses) made it clear that the system was 
not just a very bright set of high-redshift galaxies, but 
in fact was a candidate strong-lensing system. 

A literature search for known lens/source systems in 
the COSMOS field showed that this s ystem had already 
been rep orted by iFaure et all (120081 ) and iFaure et&U 
(|2011l) . IFaure et al.1 (|2008l l201ljT d"esignated the sys- 
tem as "COSMOS 0050+4901" and also n oted that it 
is a cl ear quadruple-lens in the ACS images. IFaure et al.1 
(|2011l) report a spectroscopic redshift for the lens at z spec 
= 0.960. No spectroscopic redshift was determined for 
the source galaxies which are e xtremely faint in th e ob- 
served optical bands; however, IFaure et~ai1 (|2008D esti- 
mate the redshift of the source as z p hoto = 3.34. This 
does not agree well with our photometric redshift of the 
source (§ 5.1); however, the source galaxies are extremely 
faint in the observed optical and their SEDs have no no- 
table features there (see § 5.2) so it is unsurprising that 
a highly uncertain photometric redshift was estimated 
without high-quality NIR data such as that from Ultra- 
VISTA. 

The analysis presented in lFaure etlill ()2008l l"2011[ ) fo- 
cuses on the properties of the lens galaxy. In this paper 
we focus on the properties of the source galaxy which is 
currently the only-known candidate for a strongly-lensed 
compact quiescent galaxy. 

3. LENS MODEL 

In order to measure the structural parameters of the 
source galaxy, we construct a model of the lens using 
the position and brightness of each of the four images of 
the source. We assume the lens galaxy is embedded in a 
dark matter halo with an elliptical isothermal sphere pro- 
file with an additional core component. The lens model 
has six free parameters: the Einstein radius (Rh„), core 
radius (R core ), ellipticity (e), position angle (PA), and X 
and Y position of the lensed galaxy in the source plane 
(X s , Y s ). 

In order to constrain this large parameter space, we as- 
sumed a simple Sersic profile for the source galaxy with 
an effective radius (R e ) of 1.0 kpc, and a Sersic index 
(n) of 2.5. We then generated a grid of models spanning 
the six-dimensional parameter space with fine sampling 
and used ray tracing to construct mock images at ACS 
resolution of the four images of the lensed galaxy. The 
residuals from a subtraction of the ACS image and the 
mock images were then evaluated using x 2 to determine 
the best-fit lens model parameters. The images, model 
and residual for the best-fit lensing model are shown in 
Figure 2. 

With four images of the source, the model is well- 
constrained and the best- fit parameters are R_ei„ = 

3 o4"+ 015 r - n <^"+ 0A P — n n7^+ 015 pa — - 

<J.v<± _o.ll) ^core — U.O —0 5' — U-U ' ° — 010' — 

31.0l°J degrees, X s = 0.15"t°°?, and Y s = -0.11"±°;°i- 
The error bars have been calculated using 100 Monte 
Carlo simulations where the background noise in the 
image was resampled and the lensing model was re- fit. 
The lensing model implies that the magnifications for 
the source images A, B, C, and D are factors of 4.5^6' 
5.3t^g, 4.6^5, and 2.1±g;|, respectively. 

3.1. The Size and Sersic Index of the Source Galaxy 



With the parameters of the lensing model constrained 
by the positions and brightnesses of the multiple images, 
we fix those parameters and re-perform the ray tracing 
and x 2 -fitting of the ACS images. This time allowing the 
R e and n of the source galaxy as free parameters in the 



fit. This produced a best-fit of R e = 0.43 



+2.40 
-0.40 



kpc and 



n = 1.9_ 9 . Although formally fit as an ultra-compact 
galaxy with a profile intermediate between a disk and 
bulge, the uncertainties are extremely large due to the 
low S/N of the galaxy in the ACS F814W image. 

The S/N of the ground-based K s -band image is sub- 
stantially higher than the ACS image (see §4) so we at- 
tempted to obtain better constraints on the size and pro- 
file using those data. With the higher S/N data K s -band 



data we measure an R P = 0.64 



+0.08 
-0.18 



kpc and n 



-0.9- 



This is remarkably consistent with the ACS measure- 
ment despite the much larger uncertainties on the ACS- 
measured R e . 

If we use the best-fit stellar mass of the system (§ 5.2) 
then the galaxy lies a factor of ~ 4.5 below the R e vs. 
stellar mass relation for galaxies i n the local universe de- 
fined in Ivan Dokkum et all ([20081 ). If we co mpare to the 
stellar mass-size relation at z ~ 2 defined in lBruce et al.1 
(|2012f) . we find that it is consistent with the most com- 
pact galaxies in that sample. Taken together this shows 
that the source galaxy appears to be a member of the 
massive ultra-compact galaxy population seen at z ~ 2. 

Unfortunately, at ground-based resolution (or with the 
current S/N in the ACS image) we cannot make stronger 
statements about the Sersic profile of the galaxy. We 
note that while uncertain, the best-fit Sersic index, n = 
2.2+ ' 9 , suggesting that the galaxy that is intermediate 
between a disk and a bulge. T his value is also similar to 
the median n of galaxies in the Ivan Dokkum et al.l (2008) 
sample (n — 2.3); however, we note that that sample has 
a large range of Sersic indexes (n = 0.5 - 4.5). Unfortu- 
nately, better constraints on the light profile of the source 
will not be possible until higher S/N imaging with HST 
resolution is available. 

4. PHOTOMETRIC DATA 
4.1. Deblending Method 

The lens galaxy is bright and resolved at ground-based 
resolution. Examination of the PSF-matched images 
shows that its extended light distribution causes non- 
negligible contamination to the photometry of the source 
galaxies. In order to re move this contami nation we used 
the GALFIT package ([Peng et al.l I2010D to model the 
lens/source system and subtract the flux of the lens. For 
this procedure we do not use the PSF-matched images, 
which are convolved to the worst seeing image (the g + 
image, ~ 1.2" seeing). Instead we use the original im- 
ages, which have typical seeing of < 1.0" in the optical 
bands and ~ 0.8" in the NIR band. This subtraction 
procedure is illustrated in Figure 3, and is as follows. 

In the J, H, and K s bands the source galaxies A, B, C, 
and D are clearly detected and can be distinguished from 
the lens in the images (see Figures 1 and 3). For these 
bands we used GALFIT to perform a simultaneous fit to 
all five galaxies. Given the complexity of a five-object si- 
multaneous fit, the process was somewhat unstable; how- 
ever, we continuously refined the initial guess parameters 
until a good fit with null residuals for all five galaxies was 




Model 



J v. 



Residual 




Fig. 2. — Top Panels: ACS F814W postage stamps of the four images of the source galaxy. Middle Panels: The best-fit ray-traced models 
for the source galaxy. Bottom Panels: Residuals of subtracting the model for the lens. The best-fit model simultaneously reproduces the 
relative brightness and positions of the source, as well as its morphology and size. 



achieved. In these three bands the final fit values for the 
Sersic index (n; ens ), effective radius (R e ,z e ns), axis ratio 
(qiens), and position angle (PA; ens ) of the central lens 
converged to similar values. 

Once a good fit was determined, we subtracted only 
the model for the lens from the images. The model for 
the lens galaxy and the residuals after subtraction in the 
K s and H bands are shown in Figure 3. 

In the remaining 19 optical filters as well as the Ultra- 
VISTA Y-band the source galaxies are barely detected 
(see Figure 3) and attempts at a simultaneous five ob- 
ject fit using a wide range of initial guess parameters 
with GALFIT did not return reasonable values. In these 
bands we performed only a fit to the lens. This fitting 
was stable and in all bands converged to similar values 
of ni ens , R e ,iens, qiens, and PA; e „ s . For these bands we 
subtracted the GALFIT model for the central lens (the 
Y-band subtraction is shown in Figure 3). 

4.2. Optical and NIR Photometry 

With the contribution from the lens subtracted we per- 
formed aperture photometry on the residual images using 
SExtractor in dual image mode. The source galaxies are 
brightest in the K s -band so the K s -band residual image 
was used as the detection image. Photometry was per- 
formed in 1" diameter apertures on the original seeing 
images in each band. These images have different PSFs; 
however, using the PSF-matched images for photometry 
would be sub-optimal because they have had the good 
seeing data smoothed to the larger seeing disk of the 



worse seeing image. Photometry on these images requires 
larger apertures that are more susceptible to contamina- 
tion from residual flux missed in the fitting of the lens, 
as well as flux from the other source images which are 
separated by ~ 2". We correct the photometry in each 
individual band to a total flux by assuming the source 
galaxies are point sources and extrapolating the growth 
curves of bright stars in the field. Formally this is not 
completely correct because the source galaxy is resolved 
at ground-based resolution; however, as discussed in § 3, 
the R e of the source is a factor of two smaller than the 
PSF size, so a point-source provides a reasonable first- 
order approximation to the observed profile. 

The PSFs in the UltraVISTA YJHK S are remark- 
ably similar and the seeing is 0.75", 0.76", 0.79" and 
0.82" in the K s H, J, and Y bands, respectively 
(|McCracken et afll2012| ). The correction to total mag- 
nitudes in these bands ranges between a factor of 2.04 
and 2.44, hence the NIR colors do not depend strongly 
on the aperture corrections. The P SF and seeing rang e 
in the optical is larger, 0.5" - 1.2" (jCapak et al.ll2007ft . 
The aperture corrections range between a factor of 1.79 
and 3.33. This is larger; however, as we discuss below, 
the source galaxies are not detected at > 2a in the vast 
majority of the optical bands, (the exception is the deep 
COSMOS i' filter where they are detected at ~ 4-5er), so 
these bands provide only weak constraints on the SED 
of the galaxies. 

Photometry and associated errors for the source galax- 
ies as well as the lens are listed in Table 1. Calculation 
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TABLE 1 
Photometric Data 



Filter 


Lens 


Source A 


Source B 


Source C 
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Flux Density 


Flux Density 
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1 
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1 
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1 

± 


tj 


.28 


2.00 


f 

Z 


41.25 







.39 


0.53 







.61 


0.91 
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=b 
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=t 
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.14 


2.04 
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± 


1 


.81 


0.97 


i 

± 
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0.21 


_i 
± 
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± 


1 


.81 


2.94 
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2.24 







.62 


0.07 


d= 





.82 
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0.83 


0.75 


dz 
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0.01 


± 





.83 


1.96 


7A527 


4.50 
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.68 
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0.58 


± 





.86 


0.57 


± 


0.87 


1.58 


± 


0.86 


0.26 


± 





.89 


2.94 


7A827 


31.68 


± 





.97 


0.04 


± 





.90 


1.32 


± 


0.92 


1.42 


± 


0.90 


-0.03 


± 





.96 


3.23 


Y 


44.47 


± 





.48 


1.24 


± 


1 


.10 


2.44 


± 


1.13 


2.05 


± 


1.10 


0.83 


± 


1 


.17 


2.44 


J 


46.44 


± 





.35 


2.55 


± 





.78 


3.18 


± 


0.80 


2.98 


± 


0.78 


0.77 


± 





.85 


2.27 


H 


44.78 


± 





.28 


11.06 


± 





.57 


12.00 


± 


0.58 


9.88 


± 


0.57 


5.65 


± 





.61 


2.13 


K s 


36.38 


± 





.23 


10.68 


± 





.31 


10.95 


± 


0.32 


10.23 


± 


0.31 


4.72 


± 





.34 


2.04 


3.6^m 


17.58 


± 





.07 


7.45 


± 





.27 


7.38 


± 


0.32 


7.07 


± 


0.24 


2.00 


± 





.32 


5.88 


4.5/^m 


9.85 


± 





.06 


6.11 


± 





.20 


6.52 


± 


0.21 


5.75 


± 


0.18 


3.52 


± 





.23 


5.88 


5.8^m 


5.13 


± 





.19 


3.28 


± 





.12 


4.06 


± 


0.42 


3.47 


± 


0.42 


1.27 


± 





.13 


6.67 


8.0fj,m 


1.94 


± 





.11 


1.05 


± 





.21 


1.32 


± 


0.21 


0.90 


± 


0.21 


0.52 


± 





.21 


5.88 



Note. — All flux densities are listed in units of 10 19 ergs s 1 cm 2 A 1 . 

of the photometric errors includes three sources of uncer- 
tainty. We estimate the uncertainty in the background 
subtraction by calculating the variance of the background 
in 1" empty apertures in a regions around the lens/source 
system. We also estimate a contribution from photon 
noise using the gain of the detector and the original ex- 
posure times, although we note this contribution is much 
smaller than the background variations for sources this 
faint. Lastly, we perform aperture photometry on the 
GALFIT model at the location of each source. We use 
the square root of this flux as an estimate in the addi- 
tional uncertainty from subtraction of the model. These 
three sources of uncertainty are added in quadrature to 
determine the total photometric error. The aperture cor- 
rected photometry, the photometric errors, and the aper- 
ture corrections are listed in Table 1. 

4.3. IRA C Photometry 

The S-COSMOS IRAC imaging of the source is also de- 
blended using GALFIT. The source galaxies are clearly 
detected in all four IRAC channels; however, the large 
FHWM of the IRAC PSF creates stronger blends be- 
tween the lens and source than in the optical and NIR 
bands and simultaneous five-object fits were highly un- 
stable. The problem was further complicated by the fact 
that the source galaxies are almost as bright as the lens 
in the 5.8/Ltm and 8.0/xm channels. Therefore, instead 
of a simultaneous five object fit, we used the same ap- 
proach as in the optical bands and fit only the central 
lens. Given the strong blending we left the only free pa- 
rameter in the fit as the flux of the lens, and constrained 
the position, n ;ens , R ej ; eTls , qi ens , and PA; erls using the 
best-fit GALFIT parameters from the NIR fits. This 



constrained fit was much more stable and the residuals 
at the position of the lens after subtraction of the model 
were low. 

With the lens subtracted aperture photometry was also 
performed in 1.0" diameter apertures and corrected to a 
total flux using the growth curve of bright stars. The 
aperture corrections for the IRAC photometry arc sub- 
stantially larger than the optical data, by a factor of ~ 
6 (see Table 1). Although this correction is large, the 
IRAC PSF is well-determined and quite stable so the 
correction should be well-determined. Regardless, we 
note that the inclusion of the IRAC photometry does 
not change the best-fit values of the parameters derived 
from the SED fitting (§ 5.2). It does considerably re- 
duce the uncertainties in these parameters, particularly 
on the dust extinction . This is typical for high-redshift 
massive galaxies (e.g.. iMuzzin e t al. 20093), and there- 
fore we note that any uncertainties in the IRAC pho- 
tometry will not change the interpretation of the source 
galaxy's SEDs. 

5. PROPERTIES OF THE LENSED GALAXY 

5.1. Photometric Redshifts 

Photometric redshifts (z p hoto) for both the lens and 
source galaxy images are d etermined using the EAZY 
photometric redshift code ()Brammer et al.l [20081 . All 
twenty-seven photometric bands are used in the deter- 
mination of the redshift, and the values are listed in Ta- 
ble 2. The photometry for the source galaxies is derived 
from the deblended images, whereas we have derived the 
photometry for the lens from the original PSF-matched 
K s -selected catalog. The z p / loto of the lens galaxy is 0.98 
± 0.03 which agrees well with the spectroscopic redshift 
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Fig. 3. — Left Panels: Images of the lens/source system in the K s , H, and Y bands. Middle Panels: Models for the lens determined using 
GALFIT. Right Panels: Residuals of subtracting the model for the lens. Aperture photometry is performed on these residual images. 



of z = 0.960 measured bv iFaure et~afl (|2011h . The pho- 
tometric redstarts for the images A, B, C, and D are 2.41 
± 0.13, 2.37 ± 0.18, 2.36 ± 0.15, and 2.22 ± 0.22, respec- 
tively, and are all consistent within the la uncertainties. 

5.2. Stellar Mass, Age and Dust Content 

We perform SED fitting on the photometric data to 
determine the stella r population para meters using the 
FAST fitting code (jKriek et al.l 120091) . For the SED 
fitting we assume the redshift as the best-fit z p hoto 
from EAZY. We use the stellar population models from 
iBruzual fc Chariot! (12 003) as our default model, but also 
fit to the iMarastonl (l200"5l) m o dels. We assume so- 
lar metal l icity, a ICalzetti et al.1 (|200Q[ ) dust law, and a 
IChabrierl (|2003ft initial mass function. We parameter- 
ize the star formation history as a declining exponential 
with an e-folding timescale of r. Overall we fit four pa- 
rameters, the stellar mass (M s t or ), age since the onset 
of star formation (t), r, and the V-band dust attenua- 



tion (A„). Integration of the star formation history also 
returns a star formation rate (SFR) and a star forma- 
tion rate per unit stellar mass (specific star formation 
rate, SSFR). The best-fit SEDs for both the lens and the 
three brightest sources: A, B, and C are plotted along 
with the corresponding photometry in Figure 4, and the 
best-fit stellar population parameters and associated er- 
rors are listed in Table 2. 

The lens is best-fit as an extremely massive 
(Log(M star /M Q ) = 11.4918:^1), old (LogtAge/yr- 1 ) = 
and quiescent galaxy (Log(SSFR/yr _1 ) = - 



9.4 



+0.3\ 
-0.3) 



11.701SS 



The M star derived from our photometry 
agrees reasonably well with t he value of Log(Mg tar /M^) 
= 11.64 ± 0.03 measured bv lFaure et all (|20ll . 

As Figure 4 shows, the SEDs of images A, B, and C are 
all similar within the uncertainties. This is not necessar- 
ily expected, as small perturbations in the lens potential, 
and/or a non-symmetric source galaxy can results in dif- 
ferent locations of the source being more or less magnified 
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in each image fe.g. JSeitz et~aIlfl998HSharon et al.l feo^. 
We evaluate the overall best-fit stellar population param- 
eters for the source by considering those derived from 
images A, B, C, which are the most strongly-magnified 
images. We take the median value of all parameters, but 
consider the full range of uncertainty allowed by the fits. 
Source D, which is less magnified, is excluded in this cal- 
culation because the uncertainties are much larger than 
the other three images of the source. We find that the 
source is massive (Log(M star /M© = 10. 8^0'}, corrected 
for lensing), old (Log(Age/yr _1 ) = 9.0^4), and qui- 

-1\ _ 11 /1i+0^ 



-11.41+^) with moderate 



Fig. 4. — Spectral energy distributions for the lens (upper left) and the three brightest images of the source galaxy (sources A, B, C). The 
SEDs of the sources have been measured from aperture photometry on the images after subtraction of the central source using GALFIT. 
The best-fit Bruzual & Chariot model is plotted in blue for the lens, and in red for the source galaxies. The best-fit Maraston model for 
the sources is shown in green. The photometric redshifts of all three images are consistent at z ~ 2.4. The source galaxies are extremely 
red (J - K S (AB) = 2.6), and are best-fit by (relatively) old templates with a moderate amount of dust (A„ ~ 0.8). 

data are 5.5", which is larger than the separation be- 
tween both the source and lens, and the source galaxies 
themselves. Given that individual sources cannot be re- 
solved at 24/im we performed photometry in a 7" diam- 
eter aperture to get a total flux density for the complete 
lens/source system. Within this aperture we measure a 
flux of 26.3 ± 6.3 /xJy. Based on the aperture corrections 
in the MIPS data handbook we correct this to a total 
flux of 67.3 ± I6.I nJy. 

The aperture photometry implies there is a significant 
4c7 detection at the location of the lens/source system. 
Examination of the image shows evidence for a faint ob- 
ject within the aperture; however, we note it is offset 
by a few arcseconds to the southwest of the lens/source 
system. This clear offset suggests that the lens/source 
is not the correct counterpart. There are however, no 
other plausible counterparts in the optical or NIR im- 
ages near the 24/im detection, so if the lens/source is not 
the counterpart then the true counterpart must be highly 
extincted. 

We investigate the implications if the lens/source is 
the counterpart of the MIPS source. If we assume that 
100% of the detected flux comes from the source and none 
comes from the lens, then the implied tota l SFR rate de- 
rived from the logarit hmic average of the iDale fc Heloul 
(|2002D templates fsee IWuvts et all 120081: IMuzzin et all 
2010), uncorrected for lensing, is IHI27 Ma y r_1 - The 
total Mstar in all four images, uncorrected for lensing is 
log(M star /M Q ) = 11.97^^7. This implies a hard up- 
per limit for the SSFR of the source as Log(SSFR) = 



0.5\ 



escent (Log(SSFR/yr" 

dust extinction (A„ = O.8+0 

These stellar population parameters are quite typ- 
ical of the well-studied population of m assive ultra- 
compact quiescent galaxi es at z ~ 2 (e.g.. Kriek et all 
20081: Ivan Dokkum et all 120081: IMuzzin et all l2009allbl; 
van de Sande et alll2011l:ICassata et alll20 111) . This sug- 
gests that the source is a likely to be a strongly-lensed 
member of that population. The implied SSFR is ex- 
tremely low; however, it is derived from the rest-frame 
UV/optical SED. In order to test for any potential 
obscured star formation we examined the S-COSMOS 
MIPS 24/im image of the lens/source system. 



5.3. MIPS 24nm Photometry 

The S-COSMOS 24/tm imaging is deep enough to de- 
tect galaxies with ~ 100 M Q yr _1 of star formation at 
z ~ 2.4. Unfortunately, the FWHM resolution of the 



-9- 93 -o.20- This is ~ lA dex hi S her than the SSFR de- 
termined from the SED fitting (see Table 2). If 100% of 
the flux actually comes from the source, then it implies 
that there may be some optically thick star formation 
continuing within the galaxy, or that it contains an ob- 
scured AGN, or both. We also checked if the lens/source 
was detected in either the XMM X-ray and VLA ra- 
dio observations and catalogs of the COSMOS field but 
found no detection at the location of the lens/source. 

Recent papers have shown evidence for a "main se- 
quence" of star forming galaxies out to z ~ 3 (e.g. , 
Noeske et alJl2007ljElbaz et al J 120071 ; iWuvts et al.|[201lt 
Whitaker et~aT] 12012f t. We note that while the MIPS 
flux of the source galaxies may imply it is not com- 
pletely quiescent, the typical SSFR of galaxies with 
Log( Mst ar /M m ) = 10 8 at z ~ 2.5 is Log(SSFR) = -8.6 
yr _1 (Whitaker et al. 2012). This is more than an order 
of magnitude higher than the upper limits on the SSFR, 
suggesting that if the source is a star forming galaxy, it 
lies well off the star forming main sequence and may be 
en route to becoming fully quenched. 

We also cannot rule out the possibility that some or 
all of the 24/xm flux comes from the lens. The SED- 
measured SFR of the lens is only 0.6 M Q yr -1 , which 
would only account for 4 /iJy, or 6% of the observed flux. 
However, if we note that only 10 M Q yr -1 of obscured 
star formation in the lens could account for 100% of the 
observed 24/um flux. Likewise, an obscured AGN could 
account for a fraction or all of the total flux. 

6. DISCUSSION AND CONCLUSION 

Overall, the size constraints from the lens modeling 
as well as the SED modeling confirm that the strongly- 
lensed galaxy COSMOS 0050+4901 is in fact the first 
example of a strongly-lensed massive ultra-compact qui- 
escent galaxy. The high magnification from the strong 
lensing implies that with future high S/N, high- resolution 
imaging we may now have the possibility to peer into the 
cores of these compact systems. Likewise the brighten- 
ing effects should allow us to determine a high-quality 
velocity dispersion with a reasonable integration time on 
an 8m-class telescope. 

Discovering a larger representative sample of strongly- 
lensed massive ultra-compact quiescent galaxies would 
be an important step forward to obtaining a more de- 
tailed picture of their evolution. It would also allow com- 
parison studies with the samples of strongly-lensed blue 
star-forming and submillimeter galaxies that already ex- 
ist. One question that needs to be answered in order to 
find such samples is how frequently these lensed sources 
appear, and therefore how much area with high-quality 
NIR imaging needs to be searched? This is a difficult 
question to answer precisely, as it requires knowledge 
of the underlying set of lenses and sources, as well as 
an understanding of the observational selection effects. 
Here we present an order-of-magnitude estimate based 
on what is currently known about current strong lensing 
samples, the stellar mass function, and the fact that we 
have discovered one such source in the UltraVISTA field. 

If we make the optimistic assumption that searches 
for strong lenses can identify all strongly-lensed galaxies 
down to Log(M stQr /M Q ) = 9.5 (these are the lowest-mass 
gala xies that have been identified with strong lensing, 
e.g.. lRigbv et alJl20lUlWuvts et alj|2012j) . then the frac- 



tion of high-redshift strong lenses that will be massive 
galaxies at z ~ 2 can be estimated using the stellar mass 
function. If w e integrate the stellar m ass function at 2 
< z < 3 from iMarchesini et al.l ([2009D . we find that the 
ratio of galaxies with log(M star /M0) > 11.0 compared 
to those with 9.5 < log(M sta ,./M0) < 11.0 is a factor 
of 20.7. This suggests that roughly one in twenty high- 
rcdshift strong lenses will be a massive galaxy, although 
the true number recovered will depend on observational 
selection effects. Of thes e massive galaxies, only ~ 50% 
will be a quiescent (e.g.. lKriek et al.ll2008HMuzzin et al.l 
I2009al) . 

In their "high-qualit y" sample of 16 st rong lenses in 
the COSMOS field from lFaure et al.l (|2008D estimate that 
9/16 strong lenses are at z > 2, and this is consistent 
with the median redshift of the brightest strongly-len sed 
galaxies in cluster surveys (e.g., iBavliss et all 1 20 11 alibi ). 

To estimate the frequency of strongly-lensed massive 
quiescent galaxies at z > 2, we use t he total sample of 
67 strong lenses in iFaure et al.l (|2008T) as an estimate of 
the total number of strongly-lensed galaxies per COS- 
MOS field (~ 1.5 deg 2 ). Using this empirically-derived 
number as a starting point does account for some of 
the observational selection effects. The total number of 
strongly-lensed massive quiescent galaxies that are ex- 
pected within the COSMOS field should be (total lenses) 
x (fraction of lenses at z > 2) x (fraction of lensed 
sources that are log(M star /M©) > 11.0) x (fraction of 
log(M sfQr /M ) > 11.0 galaxies that are quiescent). Fill- 
ing in the numbers above gives 67 x 0.56 x 0.05 x 0.5 
= 0.94. The number is of order unity suggesting there 
should be approximately one such source per COSMOS 
field, consistent with the single source we have discov- 
ered. This agreement is comforting, and suggests that 
the order-of-magnitude calculation is reasonable; how- 
ever, all of the numbers used in the estimate are highly 
uncertain, and most likely there are additional selection 
effects that are not accounted for. 

The true frequency of strongly-lensed massive ultra- 
compact quiescent galaxies may be factors of several 
larger or smaller. Still, the fact that we found only one 
source in the COSMOS/UltraVISTA survey would seem 
rule out the possibility that these sources are significantly 
more abundant than our estimate. Also, it is worth not- 
ing that no other such sources have yet been reported 
in sim ilar sized sur veys such as the UPS (IWilliams et al.l 
[2009h and NMBS (| Whitaker et al.ll20Tl . which would 
again support the idea that the abundance is not much 
higher than this estimate. 

It is possible that the abundance is lower than our es- 
timate, and that it is fortuitous to find such a system 
in the COSMOS/UltraVISTA field. The lack of such 
sources in the UDS, and in the 50% of the NMBS that is 
not in COSMOS field means that the current total area 
of deep- wide NIR surveys is ~ 2.5 deg 2 , with only one 
lensed massive, ultra-compact, quiescent galaxy discov- 
ered so far. It is difficult to say more without better 
statistics, but if the order-of-magnitude estimate is cor- 
rect it suggests that the space density of strongly-lensed 
massive quiescent galaxies at z > 2 is roughly one per 
1-2 deg 2 . This suggests that in order to obtain a sample 
of 10 such galaxies, medium-deep optical/NIR photom- 
etry with good angular resolution (to avoid lens/source 
blending) covering ^10-20 deg 2 would be required. Op- 
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TABLE 2 
Stellar Population Parameters 



Parameter 


Lens 


Source A 


Source B 


Source C 


Source D 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Bruzual & Chariot 2003 Models 


z photo 


0.98 ± 0.03 


2.41 ± 0.13 


2.37 ± 0.18 


2.36 ± 0.15 


2.22 ± 0.22 


Log(Stellar Mass) (M/M ) 
Log(r) (yr) 
Log(Agc) (yr) 

Log(SFR) (M yr- 1 ) 
Log(SSFR) (yr- 1 ) 


11.491°:" 

8.5t°;i 
9.4+°j 
6+ ' 3 

U.D_ 4 

-0.2ll° 9 ^ 
-11.70+°$ 


10.82+°-°f (11.48) 

-l-26± 9 , g 8 9 (-0.61) 
-12.10±°9 8 | 


10.78+8^(11.50) 

8.i'+;-; 

9 0+°' 4 
y.u_ 2 

1.0 + o 7 

-0.54l°g 1 |(0.18) 
-ll-32l° 9 ^ 


10.70t°; I (11.36) 

o 9 +0.3 
8 "'-0.3 
q 1+0- 4 

6+°' 5 
U.D_ 4 

-0.7llg;g?(-0.05) 

11 41+069 


10.68+Q "(11.00) 
7-0+o:o 

g-o+ii 

o.8l° ;I 

-99.9+. 9 5 (-99.9) 

qq q+87.60 
-yy.y_ 


Maraston 2005 Models 


Zphoto 


0.98 ± 0.03 


2.41 ± 0.13 


2.37 ± 0.18 


2.36 ± 0.15 


2.22 ± 0.22 


Log(Stellar Mass) (M/M ) 
Log(r) (yr) 
Log(Age) (yr) 

Log(SFR) (Mq yr- 1 ) 
Log(SSFR) (yr- 1 ) 


11 40+ - 00 

o 7 +0.0 
q 6 +0.0 

0+°'° 

co+o.oo 

-11 94+ ' 00 
11,3 -0 .00 


10.76+^ (11.41) 

s-2±S:l 

2+ ' 2 
-0.97+°, 7 04 (-0.32) 
12 S2+ ' 77 

±Z.OZ_gg g 


10.69l8i!(H-41) 

S.3+?; 2 

9.2±°t 

3+°' 2 
-1.64+8 9 8 J(-0.92) 

11 er+0.10 

il.OO_gg g 


10.68+°-£| (11.34) 

o q+0.2 
».3_o 2 

9 2+ 04 
a - -0.2 

2+ 02 
-0.79+° g ° 9 '(-0.13) 

11 55+0- 10 
ii.oo_ 10 


10.72+°- 1 ®(11.04) 
7 9+ 07 
q 3+0-2 

2+ - 4 
u - z -0.2 

-8.22+™ 9 (-7.56) 
-IS-S^g 



Note. — The values listed for Log(Stellar Mass) and Log(SFR) are corrected based on the lensing model. The values in brackets are 
the direct measurements with no correction for lensing. 



tical/NIR surveys with these requirements are currently 
being performed, which suggests that if searches are care- 
ful, the prospects are good for obtaining a real sample of 
strongly-lensed massive ultra-compact quiescent galaxies 
within the next few years. 
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